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Cavity cooling via quantum back-action force can extract thermal fluctuations from a mechanical 
resonator to reach the quantum ground state. Surface or bulk two-level system (TLS) defects in a 
mechanical resonator couple to the mechanical mode via deformation potential and can affect the 
cooling process significantly. Here, we develop a theory to study the cavity cooling of a mechanical 
mode in the presence of a TLS defect using the adiabatic elimination technique. Our result shows 
that the cooling process depends strongly on the resonance and the damping rate of the TLS. 



PACS numbers: 



I. INTRODUCTION 



Mechanical resonators in the quantum limit can be 
used to explore macroscopic quantum effects and the 
quantum-to-classical boundaries in such systems jA To 
reach the quantum ground state of the mechanical sys- 
tems, many approaches have been explored, including 
feedback cooling, dynamic back-action cooling, and cool- 
ing via quantum bits. 2 4 Mechanical resonators with a 
wide range of frequencies and high qualify-factors have 
been fabricated 5,6 and the interface between mechanical 
modes and other quantum systems has also been studied 
for quantum information processing^ 

Among the various approaches to reaching the quantum 
ground state of a mechanical resonator, sideband cooling 
can be achieved by coupling the resonator to an optical 
or microwave cavityA"— With the cavity driven at the 
detuning A5, the cooling (heating) rate T_ (r+) can be 
written as 



±A 6 ) 2 ] 



(1) 



where uj m is the mechanical frequency, ^0 is the cavity 
damping rate, and go is the effective linear coupling be- 
tween the mechanical mode and the cavity mode under 
the driving. The cooling rate reaches a maximum of 
T_ = 4^q//^o when the cavity detuning is at the first 
red-sideband frequency with — A5 = uo m . The cooling 
process corresponds to an energy up-conversion of the 
thermal phonons in the mechanical resonator to the cav- 
ity photons. Given the mechanical damping rate 7 m and 
the thermal phonon number n^, the steady state phonon 
number under the cavity cooling is 



(r + + 7m n^)/[(r_-r + ) + 7m ], 



(2) 



which is ultimately limited by the heating rate T+. In 
recent experiments, the resolved-sideband regime with 
<^ uom has been demonstrated which shows that it 
is promising to reach the quantum ground state in such 
systems. 
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Molecular adsorbates and crystal defects exist on the sur- 
face or in the bulk of mechanical resonators^"— These 



structures can be modeled as two-level systems (TLS) 
and their acoustic and thermodynamic properties have 
been extensively studied in amorphous solids j 22 i 23 Re- 
cently, it was experimentally demonstrated that both the 
mechanical resonance and mechanical damping can be af- 
fected strongly by the TLS defects* 2 ^! The vibration of 
a mechanical resonator modulates the asymmetric energy 
of a TLS via the deformation potential, and hence gen- 
erates a coupling between the TLS and the mechanical 
mode. In a recent work, 26 this coupling and its effect on 
the decoherence of the mechanical resonator have been 
thoroughly studied. In this work, we study the cavity 
cooling of a mechanical resonator in the presence of a TLS 
using the adiabatic elimination technique that is widely 
exploited in quantum optics. The energy spectrum of the 
coupled resonator-TLS system contains polariton states 
that are quite different from the harmonic oscillator spec- 
trum of a bare mechanical mode^ The cooling process, 
strongly depending on the energy spectrum, can hence 
be affected significantly. As we will show, the "simple" 
approach of adding the TLS and its coupling directly to 
the cooling equation of the mechanical resonator does 
not describe the cooling accurately. We derive the mas- 
ter equation for the resonator-TLS system and study how 
the steady state will be affected by the properties of the 
TLS. This theory can be extended to study the cool- 
ing of a mechanical resonator coupled with more com- 
plicated structures such as multiple TLS's. The theory 
can also be extended to include the dynamics of the TLS 
in the cooling process. Our results can be used to ana- 
lyze sideband cooling schemes for mechanical resonators 
when taking into account the defects. This paper is orga- 
nized as the following. In Sec. [Ill we present the Hamil- 
tonian of the coupled system including the mechanical 
mode, the TLS, the cavity used in the cooling process, 
and the bath modes of the thermal reservoirs for each of 
the sub-systems. Then, in Sec. HIH we derive the cooling 
master equation for the coupled resonator-TLS system 
using the adiabatic elimination technique. The results 
from this master equation and the consequence of the 
presence of the TLS on the cooling will be discussed in 
detail in Sec. lIVI Finally, discussions and conclusions will 
be presented in Sec. [Vj 
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II. THE COUPLED SYSTEM 

Our system consists of the mechanical resonator, the TLS 
defect, and a cavity that couples with the resonator via, 
e.g. optomechanical force. 2 As is shown in Fig. [TJ the 
TLS couples with the resonator with the Hamiltonian 

+ h[A z /2 + \(a + a^}a z + (hA x /2)cj x (3) 

where a (aJ) is the annihilation (creation) operator for 
the mechanical mode, a x ^ z are the Pauli operators, A z is 
the asymmetric energy of the TLS, and is the tun- 
neling matrix element between the two sites of the TLS. 
The mechanical displacement generates a strain tensor 
in the location of the defect, and hence generates a cou- 
pling between the resonator and the TLS. The coupling 
amplitude A is proportional to the deformation poten- 
tial and the second order derivative of the mechanical 
displacement and has been studied in detail in previous 
work i 22 i 26 For convenience of discussion, we rewrite the 
Hamiltonian as 

H T = huj m a^a + (hu z /2)a z + h\{aa+ + a)cr-) (4) 

in terms of the rotated Pauli matrices a XiZ with 

a z = (A z /u z )a z + (A x /u z )a x , (5) 

where uo z = yA 2 + A 2 is the frequency of the TLS and 
A = \{A x /uj z ) is the coupling constant projected in the 
rotated basis. Here, we have applied the rotating wave 
approximation to omit the term X(A z /u z )(a + o))a z and 
the counter rotating terms which only have higher or- 
der effects on the cooling process. The Hamiltonian in 
Eq. (|4]) has the form of the Jaynes-Cummings Mode l 27 i 28 
that has an energy spectrum distinctly different from the 
phonon spectrum of a bare mechanical mode as is illus- 
trated in Fig. [H When adding the cavity mode into this 
system, we consider a strong red-detuned driving applied 
to the cavity which generates an effective linear coupling 
between the mechanical mode and the cavity mode&i 2 - 
The total Hamiltonian can be written as 

H t = H T + (-HA b )tfb + Hg (a + a^)(b + 6+) (6) 

where A5 is the cavity detuning and b Qfi) is the annihi- 
lation (creation) operator of the cavity mode. 

In addition to the quantum components, each sub- system 
is subject to the environmental noise from a thermal bath 
which plays an essential role in the cooling process. For 
example, the cavity couples with a continuous spectrum 
of bath modes {bk} with the coupling Hamiltonian 

H cn = ^{c k b% + c* k blb) (7) 

k 

where c& is the coupling coefficient. The coupling with 
the bath modes induces cavity damping with the damp- 
ing rate ^o- The mechanical resonator and the TLS also 




(b) (c) 

Figure 1: (Color online) (a) Mechanical resonator couples 
with TLS defects, (b) Double-well potential model for the 
TLS. (c) Energy spectrum of the coupled system. The eigen- 
states are the polariton doublets labeled as |n±). The solid 
(dashed) arrows indicate transitions between states of identi- 
cal (opposite) polarizations. 

couple with bath modes in similar forms as that of H cn . 
At the moment, no definite theory or experimental results 
are available to accurately describe the TLS decoherence. 
For simplicity, we will model the environment of the TLS 
as a thermal reservoir with the damping rate j T . Simi- 
larly, the intrinsic damping rate of the mechanical mode 
is assumed to be j m . 



III. CAVITY COOLING OF THE COUPLED 
SYSTEM 

In this section, we study the cooling of the resonator-TLS 
system via a cavity mode driven by red-detuned source. 
When the coupling between the resonator and the TLS 
is strong, the effect of the coupling on the eigenenergy 
spectrum of the coupled system can strongly affect the 
cooling process. We will study the cooling process us- 
ing a master equation approach and derive the cooling 
equation by the adiabatic elimination technique. 



A. Eigenbasis 

The nonzero coupling between the mechanical resonator 
and the TLS modifies the eigenenergy spectrum of the 
resonator-TLS system, as is shown in Fig. QJ. Let | t, \) 
be the eigenstates of the a z operator of the TLS and \n) 
be the Fock states of the mechanical mode. The eigen- 
states of the coupled system include the ground state 
|0 i) and the polariton doublets 2 ^ 

\na)=c n a \nl) + s n a \(n-l)V (8) 

with n > 1 and a = =b. The coefficients of the eigenstates 
are 

cl = -s 7 ! = cos((5 n /2) 
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sl = c n _= sm(5 n /2) (9) 

where cos(£ n /2) = \J (uj tn + Suj)/2uj t n, Su = LO m - uj z is 
the off-resonance between the resonator and the TLS, and 
uo tn = v (5cj 2 + 4A 2 n. The eigenenergies of the polariton 
states are 

CJ na = num + (acjtn - 5u)/2 (10) 
for the n-th doublet. 

In the following, we write the quantum operators of each 
sub-system in terms of the eigenbasis of the coupled 
resonator-TLS system. For example, the mechanical an- 
nihilation operator a, which appears in the coupling be- 
tween the cavity and the mechanical resonator in Eq. (0 , 

can be written in the eigenbasis as a = with 
the operator 

6f = \(n-l)/3)(na\ (11) 

defined with the polariton states, and the matrix ele- 
ments 

4a = V^" 1 ^ + ^ls n f l sl. (12) 

Similarly, the spin operators for the TLS can be written 

aQ _ A^Aafi - f u _(n) _ n-l n 

as a± — cr^ witn a^ a — s a . 

The total master equation for this system including the 
cavity mode can be written as 29 

7 • -pna(3 

i = h [Ht > p] + Y mp+ E \^ ae P (13) 

n,oc,f3 

where C{o)p is defined as the Lindblad form for the op- 
erator o with 

C{p)p = 2op(J — po^o — o^op, (14) 

and the term 

c n a p = {n n«p + m of) + n n t fC{6^) (15) 

describes the damping between the states \{n — l)/3) and 
\na) generated by the intrinsic noise reservoirs for the 
mechanical resonator and for the TLS. Given a flat spec- 
trum for the intrinsic noise reservoirs with the mechanical 
damping rate j m and the TLS damping rate j r respec- 
tively, we can derive 

rr = l4al a 7m + k$lV (16) 

where = (exp (fujj na p /kBT) — 1) _1 is the thermal 

occupation number for the energy separation 

~ Una ~ ^(n-l)/3 (17) 

between the states \(n — l)f3) and \na). 



B. Adiabatic Elimination and Master Equation 

Let p T = Tr c (p) be the reduced density matrix for the 
coupled resonator-TLS system by tracing over the cav- 
ity mode. The cooling master equation for the reduced 
density matrix can be derived by applying the adia- 
batic elimination technique 30 to Eq. ([T3]) in the limit 
of Ko ^> go?7m 7 7T- Under strong cavity damping, only 
density matrix components with low photon numbers: 
pijnn) _ ^777 /c |p|n c ) for m c , n c = 0,1, need to be consid- 
ered. After integrating over the time variable t for the 
duration t > it can be shown that /?r 01 \ pr°\ and 

Pr 11 ^ adiabatically follow the time evolution of p^°\ The 
terms p± 1 ^ and p^ ^ are of the first order of the small 
ratio go/ 1^0 and the term p^ 1 ^ is of the second order of 
go/fro- Hence, keeping all the terms to the second order 
of go/fto, we can use the time evolution of p^ to ap- 
proximate the cooling master equation for p T . We derive 

^ = -i[H T ,p T ]/h+ £ l -\-t%*Pr (18) 

n,a,(3 

n,a,fi 

as the master equation for the coupled system. The first 
term in the equation includes the polariton Hamiltonian 
H T = a Kj na \na)(na\ with the modified polariton 
frequencies u na . The modified frequencies include small 
corrections of the second order of go/fto- The second 
term describes the decoherence due to the intrinsic noise 
reservoirs for the mechanical resonator and the TLS, as 
is given in Eq. (fT3|) . The last term in the above equation 
describes the cavity cooling in the eigenbasis with the 
cooling (heating) rates 

r;,^ = 5o«o/[^/4 + (oj na p ± A b ) 2 ], (19) 

where uo na (3 is defined in Eq. ([TT)) . 

The main differences between the above master equation 
and the standard cooling equation 8,9 are: 1. the cool- 
ing rates in Eq. ([T9|) depend on the energy differences 
Una(3 which are state-dependent and are modified by the 
finite coupling constant A, compared with Eq. (pQ); 2. the 
matrix elements which is the projection of the an- 
nihilation operator a in the polariton basis can be quite 
different from the factor \fn in the cooling equation for a 
bare mechanical mode. As we will see below, the matrix 
elements can play an essential role in the cooling 
process. 

For A = 0, it can be shown that Eq. (fT8j) recovers the 
form of the cooling equation for a bare resonator. 8,9 For 
a finite A, the cooling process can be significantly altered. 
In the dispersive regime with \(j m — uj z \ ^> A, A^2 ~ Vn 
for the transitions between states with identical polar- 
ization, but A™ = (X/Slj) and A™ = 0[(X/Suj) 3 ] ~ 
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for the transitions between states with opposite polar- 
ization. The mechanical cooling therefore includes two 
separate cooling ladders each involving states with iden- 
tical polarization, plus an additional cooling for the TLS 
with a much smaller cooling rate ~ T-(X/5uj) 2 , as is 
shown in Fig. [Tfc. This is further confirmed by an analyt- 
ical study of the cooling process in the dispersive regime 
where the cooling can be studied by applying the unitary 
transformation^ 

U = exp [-(aa+ - a f cr_)A/M (20) 

to Eq. ([T3]) . After the transformation, the TLS becomes 
decoupled from the resonator but obtains an extra cou- 
pling with the cavity mode: 

#r,c = g (\/6u)((T+a + cr_a f ) (21) 

which is to the first order of the factor X/Suj. This cou- 
pling generates cooling (polarization) of the TLS with 
a cooling rate ~ T_(X/Suj) 2 . The cooling of the me- 
chanical mode also recovers the results in Eqs. (pQ) and 
([2]) . Detailed comparison shows that this analytical result 
agrees with the numerical solution for the steady state of 
Eq. ([T8|) . In contrast, in the near-resonance regime with 
\oJm ~ oj z \ ~ 0, A^l a = (y/E ± V^l)/2 for n > 2 and 

A^2i = l/y/2. For large n, A^l a l/Ay/n, indicat- 
ing the vanishing of the transitions between states with 
opposite polarization. While for small n, the transition 
matrix elements between states with identical and op- 
posite polarizations are comparable, indicating a strong 
mixing of all the low- lying states in the eigenbasis. 



IV. RESULTS 



To illustrate the effect of the TLS on the cooling pro- 
cess, we numerically solve the steady state of Eq. ([T8]) . 
The dependences of the steady state phonon number of 
the mechanical mode and the TLS polarization on the 
TLS frequency, the TLS damping rate, and the cavity 
detuning are studied. 

The parameters we choose are comparable with what 
have been achieved in recent experiments using mi- 
crowave or optical cavities : 14 i 27 uj m = 200 MHz, mechan- 
ical damping rate 7 m /co> m = 10 -6 (Q m — 10 6 ), frequency 
of the TLS in the range of uj z /uj m G (0.5, 1.5), TLS 
damping rate in the range of jr/^m G 5 x (10 -8 , 10 -4 ), 
X/oJm = 0.05, go/(jO m = 0.05, Ko/oJm = 0.15, and cav- 
ity detuning in the range of —Ab/uj m G (0.5, 1.5). The 
coupling constants go and A are chosen by estimating 
the geometry and material properties of the resonator. 26 
Given a bath temperature of ~ 100 mK, the initial ther- 
mal phonon number can be ~ 10. 




Figure 2: (a) n ss and (b) (a z ) ss versus the ratio u z /ujm for 
7r /cj m = 2.5 x (10~ 4 , 10~ 5 , 10~ 6 ) from top to bottom. Dash- 
dotted curves: results for zero coupling A = 0; dashed curves: 
results from the "simple" approach (see text in Sec. |VJ) for 

7r — 2.5 CJ m X 10~ 4 . 



A. Dependence on TLS Frequency and Damping 

We first study the steady state phonon number n ss and 
spin polarization {a z ) ss as a function of the TLS fre- 
quency uo z . The cavity detuning is set to be at the first 
red sideband frequency with — A5 = uj m . In Fig. [2^, 
the steady state phonon number n ss is plotted. It can 
be seen that the mechanical cooling can be significantly 
degraded as the TLS frequency approaches the mechan- 
ical frequency. With a moderate TLS damping rate of 
7r/^m = 2.5 x 10 -4 and a cryogenic temperature of 
fc#T = 10 fiuJrn, n ss can be raised by nearly 50 times due 
to the presence of the TLS when uo z = u m ; while n ss is 
only slightly increased when u z = 0.5 uj m . When the two 
frequencies are in resonance, the thermal noise from the 
reservoir of the TLS can be effectively transferred to the 
mechanical resonator and causes severe heating of the 
mechanical mode. When there is a finite off-resonance 
between the two frequencies, the off-resonance generates 
an energy barrier that prevents the transfer (exchange) 
of energy between the two sub-systems. Hence, in the 
dispersive regime, the mechanical cooling is only slightly 
degraded. 

Meanwhile, the cavity cooling process extracts and 
damps the thermal noise in the TLS via its coupling with 
the mechanical mode. In Fig. [2]o, the TLS polarization 
(d- z ) ss is plotted. The TLS is maximally polarized to 
approach the state | \) when uo z ~ oj m . With the above 
parameters, we have (cr z } ss w —0.9 at 7 r = 2.5 uj m x 10 -4 . 
A finite off-resonance prevents effective extraction of the 
thermal noise in the TLS in the dispersive regime. In 
Fig. [2ji, we observe dips in the steady state phonon num- 
ber n ss when the TLS frequency approaches the mechan- 
ical frequency uj m . The appearance of these dips reflects 
the partial reduction of the thermal noise transferred 
from the TLS to the mechanical resonator when the TLS 
is highly polarized. As a result, the heating of the me- 
chanical mode by the presence of the TLS, although very 
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Figure 3: (a) n ss and (b) (a z ) ss versus the ratio j T /wm for 
oj z /(jj m = 1.3, 1,0.7 (diamond- dashed, solid, circle-solid). 



effective as the two sub-systems are in resonance, is par- 
tially reduced. 

We also study the dependence of the steady state be- 
havior on the damping of the TLS. In Fig. [3j the steady 
state phonon number n ss and spin polarization {a z ) ss 
are plotted as a function of 7 r /u; m . It can be seen that 
n ss increases sharply with the damping rate 7 r when the 
TLS frequency is in resonance with the mechanical fre- 
quency, in contrast to the much slower increase when the 
two sub-systems are off-resonance. At a damping rate 
7r/^m > IO -6 , the TLS is only strongly polarized when 
it is in resonance with the mechanical mode (the solid 
curve in Fig. [3)3). 



B. Dependence on Cavity Detuning 
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Figure 4: (a) A^/o; m and (b) n ss at cavity detuning 
versus uj z /uj m for j T /wm = 2.5 x (10 -4 , 10 -5 , 10 -6 ) (circle, 
square, triangle) respectively. Thin straight line: results from 
the "simple" approach; dashed curve: results for zero coupling 
A = 0. 



optimal value at the red sideband frequency, so that we 
have — = cj m . However, when uo z ^ uj mi the op- 
timal cooling can be reached at a cavity detuning that 



balances the two factors, i.e. 



Am) 



drifts to an interme- 



diate value between o; m and uo z . In the deep dispersive 
regime, the noise from the reservoir of the TLS is largely 
screened from affecting the mechanical resonator. The 
cooling is again dominated by the first factor so that 
— A^ m ^ — >• bj m returns to the mechanical resonance. In 
Fig. HJd, the phonon number at the optimal detuning is 
plotted as a function of uo Zl which can be compared with 
Fig. Hi- 



Next, we study the dependence of the steady state prop- 
erties on the cavity detuning. In the cavity cooling of a 
bare mechanical resonator, the optimal cavity detuning 
A^ to achieve best cooling (lowest achievable n ss ) is at 

the red sideband frequency: — A^ = uj m . In our sys- 
tem due to the presence of the TLS, this behavior can 
be altered. Using Eq. ([T8|) and studying the steady state 
by varying the cavity detuning, we numerically obtain 
the optimal detuning A^ as a function of the TLS fre- 
quency uo Zl which is plotted in Fig.H^i. It can be seen that 
the optimal detuning is shifted away from the mechanical 
resonance as the TLS frequency varies, and demonstrates 
a non-monotonic dependence on the TLS frequency. This 
dependence can be explained as a combined effect due to 
two factors. The first factor is the cooling of the mechan- 
ical resonator. Without the TLS, the cooling is optimal 
at the detuning — A5 = uj m . The second factor is the 
polarization of the TLS which is optimal when the cav- 
ity detuning is at the TLS resonance with — A5 = uj z . 
Note that the polarization of the TLS reduces the heat- 
ing of the mechanical mode by the thermal noise from 
the TLS and hence improves the cooling of the mechan- 
ical mode. When uo z = o; m , both factors reach their 



V. DISCUSSIONS AND CONCLUSIONS 

One may ask whether the above theory gives different 
results from the "simple approach" in which the TLS 
and its coupling are directly added into the cooling mas- 
ter equation for a bare mechanical resonator? To see 
the differences between these two approaches, we study 
the steady state behavior using both the "simple" ap- 
proach and the above theory. In Fig. [2j the results from 
the "simple" approach are plotted as dashed curves for 
7 r = 2.5 uj m x 10 -4 . In Fig. 0^, the results from the 
"simple" approach are plotted as the thin straight line. 
We find that the dependence of the phonon number n ss 
on the TLS properties can be very different in the two 
approaches. In particular, as seen from Fig. the op- 
timal cavity detuning — A^ in the "simple" approach 
always occurs at the first red sideband frequency; while 
in the above theory, — A^ m ^ can shift away from the red 
sideband frequency by as much as 10% of the mechanical 
resonance. Note that the results from our theory agree 
well with the results from a brutal-force solution of the 
total master equation in Eq. ([T3|) , which further confirms 
the validity of this theory. 
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We have discussed the effect of a single TLS on the cavity 
cooling and showed that the mechanical cooling can be 
degraded strongly when the frequency of the TLS falls 
within a narrow range near the mechanical resonance. 
Given the miniature size of the resonators, it is a rea- 
sonable scenario to assume that very few (e.g. one or 
two) TLS's exist in such frequency regime in the entire 
sample i 22 i 26 Furthermore, the theory developed above 
can also be extended to study the cooling of a resonator 
coupling with multiple TLS's by deriving a master equa- 
tion for the total coupled system in their eigenbasis. We 
can also extend the theory to include the dynamics of the 
TLS defects. 

In conclusion, we have studied a theory for the cavity 
cooling of a mechanical resonator coupling with a TLS 
defect. The defect, subject to thermal noise, can add ex- 
tra heating to the mechanical mode and affect the cooling 



process. We use the adiabatic elimination technique in 
the eigenbasis of the coupled resonator-TLS system and 
derive the cooling master equation for the coupled sys- 
tem. Our results showed that the cooling of the resonator 
can be significantly affected by the thermal noise of the 
TLS when the TLS frequency approaches the mechanical 
frequency. This theory can also be extended to describe 
the cooling of mechanical resonators coupling with more 
complicated quantum structures such as multiple TLS's. 
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